Table 2. A sample from the hybrid line list, the strongest 34.4 million of which are available, in electronic form, from either the CDS 
archive http://cdsweb.u-strasbg.fr/cgi-bin/qcat7/MNRAS/ or from our website http://www.tampa.phys.ucl.ac.uk/ftp/astrodata . 
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ABSTRACT 

We build an accurate database of 5200 HCN and HNC rotation-vibration en- 
ergy levels, determi ned from existing laboratory data. 20 000 energy levels in the 
iHarris e~ 3] <i2nn2bl) linelist are assigned app roximate quantum n umbers. These as¬ 
signments, lab determined energy levels and IHarris et al ] (l2nn2bD energy levels are 
incorporated in to a new energy level list. A new linelist is presented, in which fre¬ 
quencies are computed using the lab determined energy levels where available, and 
the ab initio energy levels otherwise. 

The new linelist is then used to compute new model atmospheres and synthetic 
spectra for the carbon star WZ Cas. This results in better fit to the spectrum of 
WZ Cas in which the absorption feature at 3.56 pm. is reproduced to a higher degree 
of accuracy than has previously been possible. We improve the reproduction of HCN 
absorption features by reducing the abundance of Si to [Si/H] = -0.5 dex, however, 
the strengths of the Au = 2 CS band heads are over-predicted. 

Key words: molecular data, Stars:AGB, stars: carbon, stars: atmospheres, infrared: 
stars. 


1 INTRODUCTION. 

Cool AGB carbon stars often show a strong absorp- 
tion feature at 3 micron s. This feature was discovered 
i.Iohnson fc Mended Il97(i) and later ide ntified as the C- 
H st r etch mode of HCN and/or C 2 H 2 iFav fc R.idgewa^ 
Il97fil : iRidgewav. Carbon fc Halil Il978ll . Subsequent to the 
discovery of HCN in observed carbon star spectra, it was 
found that the line opacity of HCN has a strong effect 
upon the structure of model atmospheres (lErikssorL_et_aI 
IlhS^ J^rgensen et alJll9^ . IHarris. Polvanskv fc Tennvsor 
i2002Hl . referred to as HPT, have made publicly available 
an extensive and accurate ab initio HCN and HNC linelist. 
This linelist has been used to compute model atmospheres 
and synthetic spectra that reproduce the HCN features in 
the observed spectra of the carbon stars TX Psc and WZ Cas 
and has als o allowed the ident ification off HNC absorption 
at 2.9 ^m llHarris et alJl^O.^fl . However, the HPT linelist 
has inherent inaccuracies which result in the line frequen¬ 
cies deviating from laboratory measurements by 3 cm“^ or 
more. This error is sufficiently large to be noticeable at the 
resolving power of the ISO SWS spectrometer. In particular 
IHarris et all ll2003^ fail to accurately reproduce the observed 
HCN Q-branch absorption feature at 3.56-3.62 ^m. 

In this work we assign approximate quantum numbers 


to 20 000 of the HPT ab initio energy levels. The quality of 
the HPT linelist is improved by substituting energy levels 
derived from laboratory measurements of line frequencies, 
for ab initio energy levels. This improves the frequencies of 
the lines for which experimental data is available, to spec¬ 
troscopic precision (0.1 cm“^) and allows the calculation of 
more accurate C-star synthetic spectra. These new synthetic 
spectra significantly improve the reproduction of the HCN 
Q-branch absorption feature at 3.56-3.62 ^m. 


2 ENERGY LEVEL ASSIGNMENTS 

The linelist of IHarris et al 1 J2002bll (HPT) is calculated us¬ 
ing the only exact quantum numbers for a hetronuclear tri- 
atomic, angular momentum (J) and parity. Conversely lab¬ 
oratory data is generally presented in terms of the approxi¬ 
mate vibrational quantum numbers C-H stretch (ui), bend 
(ug), C-N stretch (ua), vibrational angular momentum {1) 
and isomer (HCN or HNC). Where I can take values from 
0 or 1, if W 2 is even or odd, and up to Z ^ V 2 and Z ^ J, in 
steps of 2. Thus for a U 2 = 4 with J ^ 4, I can take values 
of 0,2,4, but for J = 2 or J = 3, Z can only take values of 
0 or 2. Figure 0 is a schematic showing the distribution of 
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even parity energy levels as a function of J and energy. New 
energy levels appear with each higher J, these are the I = J 
states. In order to match the lab determined energy levels 
with their corresponding ab initio energy levels it is neces¬ 
sary to assign the approximate quantum numbers to the ab 

initio energy level list. _ 

_ By comparison with the assignments of iBowman et afl 

the non-rotating (J = 0) vibrational states were as¬ 
signed by hand up to 10 000 cm“^ above the zero point 
energy. We then least squares fitted the vibrational expan¬ 
sion: 

3 

E{v,l)+Eo = uii(vi + di/2) 

i=l 
3 3 

+ EE Xij{vi + di/2){vj -\- dj/2) 

i=l j = l 

-I- xH (1) 

to the assigned energy levels. Where Eq is the zero point 
energy, ui are the harmonic constants, Xij are anharmonic 
constants, the degeneracies are di = da = 1 and d 2 = 2. 
HCN and HNC sta tes were fitted separat ely. For J = 0, xi 
is zero. As found bv iMaki fc Mell^ 11200111 these vibrational 
expansions give poor fits to both HCN and HNC, however 
the fits were useful in identifying and correcting mis-assigned 
states. 

A somewhat different assignment strategy is required, 
for J > 0, as the I > 0 states do not occur for J = 0. The 
distribution of energy levels in a vibrational band of HCN 
can be approximated by using the rotational expansion for 
a linear triatomic molecule: 

E{v,J) = Eiv) + B4J{J + l)+l^] 

-D4J{J + l)+l^f 

-t- 1)) -t- / -j- ... (2) 

where E(v) is the vibrational energy of the band, is the 
rotational constant, is the centrifugal distortion constant 
and is a higher order constant. To estimate the J = 1,1 = 
0 energy levels, we used the ground vibrational state values 
of the rotational constant (B«) of 1.478221 and 1.512112 
cm“^ for HCN and HNC respectively, and the first two terms 
of eq. 121 By comparing these estimated energy levels with 
the ab initio levels the J = 1, Z = 0 states were assigned. The 
remaining states were those with 1 = 1, the lowest energy 
of these states is the (0,1^,0), bending fundamental which 
was assigned. The vibrational expansion eq.^was used to fit 
the newly assigned J = 1, Z = 0 and (0,1^,0) states and thus 
estimate values for the higher J = 1, Z = 1 states. Again, by 
comparison of the estimated values with the ab initio values, 
the J = 1, Z = 1 states were assigned. 

Higher J states were assigned in the same way, the rota¬ 
tional expansion was fit to the previously assigned data for 
each band, the Z < J states assigned by comparison. These 
assigned Z < J states were fitted with the vibrational ex¬ 
pansion, allowing the assignment of the Z = J states. This 
process was carried out on a J by J basis until energy levels 
up to J = 60 had been assigned. In total we assigned ap¬ 
proximate quantum numbers to around 20 000 energy levels, 
for HCN this is up to 3 quanta of H-C stretch, 13 quanta of 
bend and 4 quanta of C-N stretch. 


The distribution of the line strengths of a given band 
can be approximated by Honl-London factors. By fitting the 
line transition dipoles of each band using Honl-London fac¬ 
tors is is possible to identify lines which have an intensity 
which is inconsistent with that of the band as a whole. In this 
way we cross checked the consistency of our assignments, us¬ 
ing line intensities, and corrected our assignments where nec¬ 
essary. However, there are states in the linelist of HPT which 
are close enough in both energy and quantum numbers for 
intensity stealing to break the distribution of inte nsities in a 
band, see lHarris. Polyansky fc TennvsonI ||2002^. For these 
states the assignments of approximate quantum numbers, 
based on energy and/or transition dipole is difficult and ar¬ 
guably meaningless, such states are left unassigned. 

3 ENERGY LEVELS DETERMINED FROM 
LABORATORY DATA. 

HCN and HNC have been studied in detail in the laboratory 
iSmith et alJlTo^ iMaki et al.l 1200(1 iLecoutre et ffi]l200(l : 

iNorHmii^l^dTl^OTnMAk^^MeU^l^nih . HCN has been 

studied up to a maximum of 6 quanta in the bending and H- 
C stretch modes and 4 quanta in the C-N stretch mode. Due 
to its nature, HNC is more difficult to produce in the lab 
thus data is less extensive than for HCN, up to 2 quanta of 
bend, 4 quanta of H-N stretch and 2 quanta on C-N stretch. 
This laboratory data is presented in varying formats, from 
line frequencies to fits of polynomial expansions in angu¬ 
lar momentum, such as eq. |21 Using the electron i c dat a 
base of laboratory lin e frequencies from iMaki e t alJ i200(il : 
iMaki fc Mellaul (1200 jl : lNorthruD et al.l il997fl we have com¬ 
puted and compiled a list of HCN/HNC energy levels rela¬ 
tive to the zero point energy of HCN. For any given band 
some or many line frequencies remain unmeasured, there are 
therefore missing energy levels in many of the vibrational 
states. These missing energy levels can be interpolated by 
means of the rotational expansion for a linear molecule eq. 
121 For a given vibrational excitation, Z, and parity we least 
squares fit eq. |21to our laboratory determined energy levels, 
to obtain E(v), R„, and JT„. These constants are then 
used to interpolate any missing energy levels. 

In the high temperatures of a stellar atmosphere the 
HCN molecule can be excited to a very high angular mo¬ 
mentum. HPT found that at 3000 K states up to J = 60 
contribute up to 93% of the rotationally converged parti¬ 
tion function. The maximum rotational excitation in our 
laboratory determined linelist is J = 55, the minimum is 
J = 20, to match HPT we must extrapolate to J = 60. As 
the rotational expansion (eq.|21l is divergant, it is not physi¬ 
cally realistic and is therfore unsuitable for extrapolation to 
high J. In order to extrapolate to J = 60, we correct the 
ab initio energy levels of HPT by adding a constant for a 
given vibrational state, Z, and parity. This constant is given 
by: C = EiabHmax'} Eai{dmax'), where EiahH'ixiax'} is the 
energy of the state in the laboratory energy level list with 
the highest angular momentum (Jmax), and Eai{Jmax), is 
the energy of the ab initio energy level with J = Jmax- 

The interpolated and extrapolated list of laboratory 
determined energy levels is referred to as the lab/empirical 
energy level list. It contains 5200 HCN/HNC energy 
levels up to a maximum of 6 Quanta in bend, 3 quanta 
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J=0 J=1 J=2 



Figure 1. A schematic showing the distribution of even parity HCN energy levels as a function of angular momentum (J) and Energy. 
States are labelled {vi ,V 2 ,l,V 3 ). 


of H-C stretch and 3 quanta of C-N stretch and to a 
maximum angular momentum (J) of 60. We have in¬ 
corporated the lab/empirical energy levels, the 20 000 
assignments in to the original ab initio energy level list of 
HPT. A sample of this new energy level list is given in 
table 0 the full list can be downloaded from our website 
(http://www.tampa.phys.ucl.ac.uk/ftp/astrodata ) or the 
CDS archive (http://cdsweb.u-strasbg.fr/cgi-bin/qcat7/MNRAS/ 
The column labelled index in tabled is the index num¬ 
ber given to the energy level by HPT, J and P are the exact 
quantum numbers of angular momentum and parity, n is the 
number of the energy level in the J-P symmetry block, E^i is 
the value of the ab initio energy level, iso labels the state as 
either HCN (iso=0) or HNC (iso=l), vi, V 2 , I, and W 3 are the 
approximate quantum numbers, Eiab is the lab/empirical en¬ 
ergy, label is a single character label which is either ’e’ for 
a laboratory determined energy level, ’c’ for an interpolated 
energy level or’t’ for a corrected ab initio energy level. The 
error on the energy level takes three forms, for a lab de¬ 
termined energy level this is the compound error of the line 
frequency measurements used to derive the energy level. For 
an interpolated energy level this is the standard deviation 
on the fit of eq. For a corrected ab initio energy level the 
error is the difference between the energy predicted by the ht 
of eq .|21 to lab and empirical energy levels, and the corrected 
’ab initio’ energy. The errors of the lab determined energy 
levels are typically of the order of 10 “"^ cm“^, the errors on 
the interpolated energy levels are typically of the order of 
10 “^ cm”*^, and of t 
mitio energy levels. 


the line frequencies are computed from the original ab initio 
energy levels. However, if lab/empirical energy levels have 
been determined for both the upper and lower states, we can 
use the lab/empirical energy levels to compute a more ac¬ 
curate line frequency. The resulting linelist significantly im¬ 
proves the accuracy of the most intense transitions between 
the low lying energy levels. This improvement in the linelist 
Is evident in hgures and which show opacity functions 
calculated with both the ab initio and hybrid linelists. The 
predominant features in these opacity functions are the Q 
branches of the {Av 2 = 1 - 1 - Avs = 1 ) and (Ai ;2 = 1 ) bands, 
respectively. As the Q branches have a large concentration 
of lines over a small frequency range, it is with these Q 
branches that the improvements within the synthetic spec¬ 
trum are most evident. 

To maintain accuracy and consistency when computing 
a line frequency it is important to use either ab initio or 
lab/empirical energy levels for both the upper and lower 
states. The error on a given ab initio energy level in general 
increases with an increasing value of a quantum number. For 
example the deviation of the (0,0°,1) J = 0 ab initio energy 
level from laboratory determination is -1-3.7 cm“^, for the 
(0,0°,2) and (0,0°,3) states the deviation is -1-8.4 and -1-13.8 
cm“^. Thus when calculating frequencies between excited 
states using the ab initio energy levels there is, more often 
than not, a partial cancellation of error. This is not the case 
if a combination of lab/empirical and ab initio energy levels 
are used to calculate frequency, which results in a larger 
error on frequency. 


4 THE NEW LINELIST 

The lab/empirical and ab initio energy levels, together with 
the Einstein A coefficients from the HPT linelist, allow the 
generation of a new hybrid linelist. For the majority of lines, 


To reduce the required storage space and increase the 
speed of computations using our linelist, we have truncated 
the hybrid linelist at a minimum intensity of 3 x 10“^® 
cm molecule”^ at 3000 K. This results in a linelist of 
34.4 million lines less than 10% of the 393 million lines in 
the origonal linelist. These 34.4 million lines account for 
more than 99.99% of the opacity of the full linelist at a 
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Table 1. A sample from the lab/empirical energy level list, which is available in full, in electronic form, 
from either the CDS archive (http://cdsweb.u-strasbg.fr/cgi-bin/qcat7/MNRAS/1 or from our website 
Ihttp://www.tampa.phys.ucl.ac.uk/ftp/astrodata I. 


Index 

J 

P 

n 

Eai (cm 1) 

iso 

VI 

V2 

i 

^^3 

Eiab (cni error (cm 

label 

1772 

3 

1 

22 

4704.012753 

0 

1 

2 

0 

0 

4702.019193 

3.00E-04 

e 

1773 

3 

1 

23 

4719.910847 

0 

1 

2 

2 

0 

4716.909245 

4.58E-04 

e 

1774 

3 

1 

24 

4880.708077 

0 

0 

7 

1 

0 




1775 

3 

1 

25 

4909.459579 

0 

0 

1 

1 

2 

4895.763446 

2.00E-04 

e 

1776 

3 

1 

26 

4909.583711 

0 

0 

4 

0 

1 

4906.162828 

2.04E-01 

c 

1777 

3 

1 

27 

4910.798531 

0 

0 

7 

3 

0 




1778 

3 

1 

28 

4925.183847 

0 

0 

4 

2 

1 

4920.596764 

8.51E-02 

c 

1779 

3 

1 

29 

5203.782210 

1 

0 

0 

0 

0 

5203.779963 

5.00E-07 

e 

of 3000 

K. 

A sample of the linelist is given 


and for SI units {v in m’ 

^ and Ai/ in s the 


in table here v is frequency, E" is lower state energy, 
Ai/ is the Einstein A coefficient, J" and J' are lower state 
and upper state angular momentum quntum numbers, 
p is parity where 1 is even and 0 odd, n is the number 
of the level in the J-parity block, index is the unique 
label of the energy level, iso labels a HCN state if 0 or a 
HNC state if 1, ui, V 2 , I, vs are the approximate quantum 
numbers. Where the approximate quantum numbers have 
not been assigned a value of -2 is given. The full ver¬ 
sion of the linelist is available from either the CDS archive 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7/MNRAS/ or from 
our website http://www.tampa.phys.ucl.ac.uk/ftp/astrodata 
The file size is 3.6 Gb uncompressed and 821 Mb compressed 
using bzip2. The linelist is sorted is ascending frequency 
order. 

The untruncated linelist can be obtained by download¬ 
ing the original HPT Einstein A coefficients, the new as¬ 
signed energy level list and running the supplied EORTRAN 
utility program dpsort-v2.0.f90. The Einstein A coefhcient 
file from HPT is sorted using ab initio frequencies, if lab de¬ 
termined frequencies are used the linelist will no longer be 
in full frequency order. 

The integrated absorption intensity of the lines can be 
computed from the Einstein A coefficient using 


, C(2J' -b 1) 

^ 

QvrV^ 


-E" 

'W 


1 — exp 


e) 


Aif 


(3) 


where I is integrated intensity, is the ro-vibrational 
partition function, J' is the upper state rotational quan¬ 
tum number, v is the transition wavenumber, E" is the 
lower state energy, C is a constant and Aif is the Ein¬ 
stein A coefficient. To return I in cm per molecule, with 
Aif in s“^ and v in cm“^ C has the value of (87rc)“^ = 


of the llBarber. Harris fc TennvsonI ll2002^ rotationally con¬ 
verged HCN/HNC partition function. Dimensionless oscil¬ 
lator strengths are often used by those involved in the mod¬ 
elling of stellar atmospheres. Weighted oscillator strengths 
{gf) can be computed from the Einstein A coefficients using: 


^ K{2J' + 1) ^ 

gf = - -Aif 


(4) 


where A is a constant which for cgs units {v in cm ^ and 
Aif in s“^) has the value 


K = = 14991.97 s m~^ 

zTre^ 


( 6 ) 


K = = 1.499197 s cm-2 


(5) 


5 COMPUTATION OF MODEL 

ATMOSPHERES AND SYNTHETIC 

SPECTRA 

We computed th e model atmosphe res of C-giants with the 
SAM 12 program [ |^eiffidj2003^|3) which is a modifica¬ 
tion of ATLAS 12 iMiruc3ri99ijl . SAM 12 is a plane parallel 
code which makes use of the LTE approximation, hydro¬ 
static equilibrium and the conservation of flux (there are no 
sources or sinks of energy within the atmosphere). The pho¬ 
tospheres of C-giants lie on the upper boundary of their con¬ 
vective envelopes, thus convection plays an important role 
in these atmospheres. We use the mixin g length theory of 
ID convection modified bv lKurnc^ ill 999^ in ATLAS12, and 
adopt a value of the mixing length parameter of l/H =1.6. 

SAM12 uses the standard set of continuum opaci¬ 
ties included in ATL AS12, with the addition of con tinu¬ 
ous opacities from C“ iMverscough fc McDoweli|ll966^ . H^ 
llDovlel 119681), and bound-free a bsorption of C, N, O from 
TOPBA SE l^e^oi^t aU ^92() with cross-sections from 
IPavlenko fc ZhrLkovsO ~ il20033b These bound-fre e opacity 
tables have been made available on the Web jPavlenkol 
l2003(j) . Atomic and molecular lines opacities w ere accounted 
for b y using the opacity sampling approach llSneden et all 
^7^. The data for these lines was obtained from several 
sources. HCN/HNC lines taken from this wor k, atomic lines 
from the VALD database l|Ru£kae^km^99|^jCN, C 2 , SiH, 
MgH, CH from CDRQM 18 o flKurucd (Il99,‘^ CO from 
iGoorvitchl (|W93), and CS from IChan^a'e^iTlll99!^^ . The 
following molecular electronic bands were accounted for: 
CaO (C^E - X^E), CS/A^E - X^E), SO (A^H - X^E), SiO 
(EiE-X^E), SiO/A^H - XiE+, NO (C2n,- X2n,), NO(B2n, 
-X2no, no(a 2 e+ - X2n,), MgO(BiE+ - x^e+z, Aio(c^n 

-X^E), A10(B2E+ - X^E+Z, in the framework of the JOLA 
approach. Eor WZ Cas these electronic bands are only week 
sources of opacity. For all cases we adopted the carbon iso¬ 
topic ratio of 5. The shape of every line was deter¬ 

mined using the Voigt function H{a,v). Damping constants 
were taken from line data bases or computed using Unsold’s 
approach iUnsold Ill95,'^^ . In model atmosphere computa¬ 
tions we adopt a microturbulent velocity, Vt, of 3.5 km/s. 

Model atmospheres were computed using the ’best fit’ 
effective temperature and C/O ratio, found in our earlier 
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Wavelength [jim] 

Figure 2. The opacity function of HCN/HNC at 2800 K, the lines have been broadened by a Gaussian with half width at half maximum 
of 0.003 /im. 



Wavelength [|im] 


Figure 3. The opacity function of HCN/HNC at 2800 K, the lines have been broadened by a Gaussian with half width at half maximum 
of 0.003 fim. 


work iHarris et alj 1200 . Teff= 2800 K, log(Nc/No) = 
—0.003, log(g) = 0.0. Figure 21 shows the pressure temper¬ 
ature curves of model atmospheres computed with the ab 
initio linelist and hybrid linelist. There is little difference be¬ 
tween the atmospheres at depth, but at lower optical depths 
the model atmosphere computed with the hybrid linelist is 
slightly hotter at a given pressure. 

S ynthetic spectra were calculated with the WITA6 pro¬ 
gram (|Pavlenkoll200nh . using the same approximations and 
opacities as SAM12, but with a considerably hner grid of 
wavelength points (AA = 0.01 — 0.02 A). We have also 
computed synthetic spectra with reduced S and Si abun¬ 
dances of [Si/H]=0 and -0.5 dex, to attempt to better ht 
HCN absorption and reduce the strength of CS absorption, 
see section O 


6 COMPARISON OF SYNTHETIC SPECTRA 
WITH OBSERVED STELLAR SPECTRA. 

The observations of WZ Cas used in this study were taken 
as part of the Japanese gu aranteed observing tim e (RED - 
STARl, PI T.Tsuji, e.g., lAoki. Tsuii fc Ohnakal lllOOSh l. 


The data and data reduction used in this work is discussed 
in detail in our earlier paper iHarris et al lliooi). 

Synthetic and observed spectra of WZ Cas between 2.7 
and 4.0 microns are shown in hgure 21 At this resolution 
there is little difference between the synthetic spectra cal¬ 
culated with the ab initio and hybrid linelist. A synthetic 
spectrum computed with solar abundances results in HCN 
absorption that is stronger than in the observed spectrum 
WZ Cas. We find that [Si/H]=-0.5 provides the best ht of 
HCN features. The reason for this is that in carbon rich 
atmospheres silicon and sulphur have a strong second ary af¬ 
fect o n chemical equilibrium, after carbon and oxygen (iTsuii 
[^T^). Indeed, SiS is analogous to CO and is a very stable 
molecule, thus atmospheres with a high ratio of Si to S are 
silicon rich and those with a low ratio are silicon poor. Figure 
21 shows the number densities of CS, SiO, SiS and HCN, as 
a function of pressure in our model atmospheres computed 
with [Si/H]=0 and -0.5 dex. As the abundance of silicon 
is reduced from [Si/H]=0 to -0.5 dex, there is less SiS and 
SiO, freeing S and O to form CO and CS, this reduces the 
abundance of HCN. 

At higher resolutions the improvements in line positions 
in the hybrid linelist become more obvious. Figures |71 and 
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Figure 4. Temperature pressure curves of model atmospheres of Teff = 2800 K, log(Nc/No) = —0.003, log(g) = 0.0, computed with 
the ab initio and hybrid linelists. 


1 ^ show the synthetic spectra calculated with the ab initio 
and hybrid linelists, and the observed spectrum of WZ Cas, 
between 3.55-3.61, and 3.84-3.895 fj,m respectively. There is 
a clear and improvement of the fit to observation, across the 
whole of the 3.55-3.61 fim range. The remaining discrepan¬ 
cies between observed and synthetic spectra over this range 
are likely to be due either to the spectra of missing chemi¬ 
cal species or to transitions between states with V 2 > 5, for 
which laboratory data is unavailable. In figure[7| we have la¬ 
beled the position of the band centres of several HCN bands. 
Absorption features can be seen close to the band centres, 
these are due to the Q-branch (AJ = 0) lines of the band, 
which are closely spaced and result in strong absorption over 
a narrow range. 


6.1 The CS An = 2 band heads and the abundance 
of Sulphur and Silicon. 

In the 3.84-3.895 fim region shown in figure|Hl we have iden¬ 
tified the Q-branches of the 3 lowest lying HCN states. The 
reproduction of the shapes of these HCN absorption features 
is significantly improved, by the use of the hybrid linelist. 
However, we over-predict the st rength of the CS v = 0 —> 2 
band head at about 3.868 /rm. lAoki et alJ ^1998^ had the 
same problem, they suggested that the poor fit was a re¬ 
sult of emission by CS from a circumstellar shell. We feel 
that this is an unlikely scenario. Absorption from both the 
fundamental and the hot bands is clear in our synthetic spec¬ 
tra, and also appear, all be it more weakly, in the observed 
spectrum. In contrast, emission from CS at temperature sig¬ 
nificantly lower than that of the stellar photosphere, will, be 
predominantly from the fundamental and the low lying hot 
bands. Thus the signature of atmospheric absorption with 
emission from a cool circumstellar shell would be a weak 
fundamental band and strong hot bands, this is not seen in 
the observed spectra. 

As an alternative to circumstellar shell emission 
wetested alternative metal mixes. The problem, however, 
cannot be solved by reducing the sulphur abundance, be¬ 
cause not only are the strength of the CS lines reduced, but 
chemical equilibrium is altered resulting in a greater abun¬ 
dance of HCN, so that the HCN features from 2.8-4 /rm 


were poorly reproduced. Similarly increasing the abundance 
of silicon to [Si/H]=0 dex weakens slightly the CS features, 
but increases the stength of HCN absorption. The net ef¬ 
fect is little improvement in the reproduction of absorption 
over the 3.8-4 micron region, see figure 1^ but a worse repro¬ 
duction of the HCN absorption accross the wider spectral 
range. 

We use the CS linelist of IChandra et all (ll995^ in our 
computations, we do not believe that errors in this linelist 
contribute to the over prediction of the strength of CS fea¬ 
tures in our synthetic spectra. This CS linelist was ca l culate d 
by using Dunham coefficients from JWinke^^l ] il984ll 
and t he dipole moment surface of iBotschwina fc Sebaldl 
1 ^8^. Within the rigid rotor approximation, the rotational 
J =0^1 transition dipole is equal to the permanent dipole 
moment, thus can be used as a rough check for the accu¬ 
racy of intens ity calculations. The J =0—>1 transition dipole 
computed bv IChandra et ahl ( Il995|) is 1.889 D, the perma- 
nent dipole moment quoted bv lBo^hwina fc Sebaldl (|^8^ 
is 1.958 D. The vibrat i onal b and transition dipoles given by 
IBotschwina fc Sebaldl ^1985^ are 0.1580 D and 0.00911 D 
for the n = 0 ^ 1 and n = 0 —> 2 transitions respec- 
tiv ely. The n = 0 —> 1 and v = 0 —> 2 transitions given 
bv IChandra et ahl (ll995^ are 0.1569 D and 0.00788 D. Al¬ 
though there are small differences, this is not large enough 
to account for the difference between observation and syn¬ 
thetic spectra. A final possibility is that there is an opacity 
source present in the stellar atmosphere which we have not 
accounted for in our models. This could be from a species 
such as C 2 H 2 , which is known to exist in some C-star atmo¬ 
spheres, but is not included in our calculations, as no C 2 H 2 
linelist is available. It is possible that increased opacity from 
such species could affect both the structure of model atmo¬ 
spheres and/or change the synthetic spectra. 


7 CONCLUSION. 

We have assigned approximate quantum numbers to 20 000 
of the ab initio e nergy l evels found in the HCN/HNC linelist 
of lHarris et al ] ll2002tJl . Existing laboratory measurements 
of HCN and HNC line frequencies are used to accurately 
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Figure 5. Synthetic and observed spectra of WZ Cas. The synthetic spectra have been convolved with by a Gaussian with half width 
at half maximum of 0.003 pm. The synthetic spectra have constants of 0.2 and 0.4 subtracted from the flux. 



Figure 6. Number densities of HCN, CS, SiO and SiS as a function of pressure in the model atmospheres with [Si/H] = 0 and -0.5 dex. 



Figure 7. Synthetic and observed spectra of WZ Cas. The synthetic spectra have been convolved with by a Gaussian with half width at 
half maximum of 0.003 pm. The synthetic spectra have constants of 0.1 and 0.2 subtracted from the flux. HCN band centres are labelled 
with arrows and an index number for the band which is tabulated in table 
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Table 3. Laboratory band centres for bands with Ai ;2 = l,Ar ?3 = 1. 


index ^ 

,v'^) 

«,u'',r,u'') 

Ac ^ 

1 

(0,1,1,1) 

(0,0,0,0) 

3.5643 

2 

(0,2,0,1) 

(0,1,1,0) 

3.5841 

3 

(0,2,2,1) 

(0,1,1,0) 

3.5652 

4 

(0,3,1,1) 

(0,2,0,0) 

3.5845 

5 

(0,3,1,1) 

(0,2,2,0) 

3.6040 

6 

(0,3,3,1) 

(0,2,2,0) 

3.5659 

7 

(0,1,1,2) 

(0,0,0,1) 

3.5953 

8 

(0,4,0,1) 

(0,3,1,0) 

3.6041 

9 

(0,4,2,1) 

(0,3,1,0) 

3.5849 


^ This is the index of the band used in figure 171 
^ From the laboratory measurements of iMaki et al.l 



Figure 8 . Synthetic and observed spectra of WZ Cas. The synthetic spectra have been convolved with by a Gaussian with half width 
at half maximum of 0.003 pm. 


determine 5200 energy levels. Both the energy level assign¬ 
ments and the lab/empirical energy levels h ave been in¬ 
cluded along with the ab initio energy levels of lHarris et alJ 
J2002tJ^ in a new energy level list. A combination of the 
lab/empirical and ab initio ene r gy leve ls are used in con¬ 
junction with the iHarris et al.l 1120021^1 Einstein A coeffi¬ 
cients to produce an improved linelist. We make the new 
energy level file which includes assignments, laboratory de¬ 
termined energy levels and the resulting improved HCN and 
HNC linelist publicly available. 

The new linelist has been incorporated into our compu¬ 
tations of C-rich model stellar atmospheres and synthetic 
spectra. The new synthetic spectra ar e a clear improve - 
ment over our earlier synthetic spectra iHarris et alJl20031 . 
and show better agreement with the observed spectrum of 
WZ Cas. The improvement is clearest in the range 3.56-3.62 
^m. We further improve the fit to HCN features by reduc¬ 
ing the Si abundance to [Si/H]=-0.5 dex, however the CS 
absorption feature at 3.8-4.0 microns is still poorly repro¬ 
duced. 
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Table 2. A sample from the hybrid line list, the strongest 34.4 million of which are available, in electronic form, from either the CDS 
archive http://cdsweb.u-strasbg.fr/cgi-bin/qcat7/MNRAS/ or from our website http://www.tampa.phys.ucl.ac.uk/ftp/astrodata . 


V (cm 

E"(cm-i) 

J" 

P" 

n" 

J' 

P' 

n' 

Ai/ (s 1) 

index" 

index' 

iso" 

Vl” 

V2' 

r 

V3" 

iso' 

Vl' 

V 2 

1' 

V3' 

label 

2.956432 

0.000000 

0 

1 

1 

1 

1 

1 

2.263E-05 

1 

99631 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

lb 

2.957251 

12055.741477 

0 

1 

92 

1 

1 

171 

1.794E-05 

92 

99801 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

ai 

2.957637 

4888.417867 

0 

1 

10 

1 

1 

17 

2.015E-05 

10 

99647 

0 

0 

4 

0 

1 

0 

0 

4 

0 

1 

lb 

2.957871 

12150.505903 

44 

0 

161 

44 

1 

206 

6.785E-07 

141511 

140356 

1 

0 

8 

6 

0 

1 

0 

0 

0 

2 

ai 

2.958350 

11744.486505 

0 

1 

87 

1 

1 

162 

2.200E-05 

87 

99792 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

ai 

2.958460 

8691.662400 

0 

1 

37 

1 

1 

68 

1.923E-05 

37 

99698 

0 

1 

8 

0 

0 

0 

1 

8 

0 

0 

ai 



